Improved search for z^^ — > z/g oscillation in a long-baseline accelerator experiment 
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We performed an improved search for v^^ — > oscillation with the KEK to Kamioka (K2K) long- 
baseline neutrino oscillation experiment, using the full data sample of 9.2 x 10^^ protons on target. 
No evidence for a Ve appearance signal was found, and we set bounds on the Ue oscillation 

parameters. At Am^ = 2.8 x 10~^ eV^, the best fit value of the K2K disappearance analysis, we 
set an upper limit of sin^ 2S^e < 0.13 at 90% confidence level. 

PACS numbers: 14.60.Pq,13.15.+g,25.30.Pt,95.55.Vj 



Introduction. — We describe a search for Vf, appear- 
ance in a beam of v^, which is the signature of a non- 
zero value of the unknown neutrino mixing parameter 
013. In the current picture of neutrino oscillations, three 
flavors of neutrinos are related to three mass states by 
the Maki-Nakagawa-Sakata matrixQ. The mixing can 
be described by two Am^ parameters (Am^^jj^, Atoq), 
three mixing angles (6123 ~ 6'atm, 6*12 ~ 6*©, and 6'i3) 
and a CP violating phase {6). Atmospheric v oscillations 
measured by several experiments |M ^ IS and con- 
flrmed by the beam experiment are well-described by 
J'u ^ oscillations with parameters sin^ 20atm > 0.92 



and 1.5 x 10 



< ^ml,,^ < 3.4 X 10- 



eV^ Solar 



t'e — > t'/j.T oscillations with parameters in the range 
0.2 < sin^'6'0 < 0.4 and 7 x 10"^ < Ato| < 9 x IQ-^ eV^ 
are also consistent with multiple observations^^, 1^ , and 
are confirmed by disappearance of reactor De^- As yet, 
very little is known about either 613 or 6, although lack 
of observed disappearance of reactor De over a few km 
baseline^S has shown that 6*13 must be smaller than 12° 
at the Am^t^j region reported by K2K0. In a 2-flavor 
approximation, the probability of appearance of v^. in a 
beam of is given by 

^ i^e) = sin^ 26i^eSin2(1.27Am2^L/£:), 

where L is the baseline in km, E is the neutrino energy 
in GeV and Arrif,„ is in eV^. For the case that Am| 



Am^t^, > Ami 



fie 
12 



Atoq , and sin^ 20atm 



23 



1, we can 



A'^atm and sin 26'pe 



sin^ 023 sin^ 26*13 



take Amj^ 
isin^ 2^13. 

Experimental apparatus and data sample. — The KEK 
to Kamioka (K2K) long-baseline neutrino experiment 
comprises a 98% pure i^^ beam with mean energy of 
1.3 GeV, created at KEK's proton synchrotron and sent 
250 km to the Super-Kamiokande (SK) detector [ll|. 
The beam is created by colliding primary protons of 
12.9 GeV/c on an aluminum target, focusing the re- 
sulting secondary pions by two electromagnetic horns, 
and letting pions decay in a decay pipe. The contam- 
ination of Ve in the beam (beam j/g) is ~ 1% at KEK 
site. Near neutrino detectors are employed at a 300 m 
baseline to measure the beam. The near detector com- 
plex includes fine-grained detectors (FGD) and a 1 kilo- 
ton water Cherenkov detector (IKT). More details of the 
subdetectors are found in references 0, 0, Q| . 



There were two running periods in K2K: K2K-I [4.8 x 
10^^ protons on target (POT)] corresponds to dates for 
which SK was instrumented with 11,146 inner detector 
(ID) photomultiplier tubes (PMTs). K2K-II (4.4 x lO^^ 
POT) corresponds to dates after SK was rebuilt with 47% 
ID PMT density following an accident. The K2K-I period 
is divided into K2K-Ia (June 1999) and K2K-Ib (Novem- 
ber 1999 to July 2001). For K2K-Ia, the horn current 
and target diameter were different from the subsequent 
phases of the experiment, and the differences of neutrino 
fluxes and systematic uncertainties are taken into ac- 
count. For K2K-II (January 2003 to November 2004), 
the lead-glass calorimeter (LG) in FGD was replaced by 
a fully-active scintillator detector (SciBar). We have an- 
alyzed the entire K2K data sample which corresponds to 
9.2 X 10^^ POT in total. The statistics are almost the 
double of our previous published results of K2K-I p^ . 
The results in this paper are obtained with a revised sig- 
nal selection and an improved sensitivity. 

Event selection. — Our search for —>■ signal is 
based on detection of charged current quasi elastic (CC- 
QE) interaction of Ve in an oxygen nucleus: v^+n — > e+p. 
Typically the momentum of a recoil proton from CC- 
QE interaction is below Cherenkov threshold in the wa- 
ter. Therefore only the electron, visible as a single-ring, 
shower-type event, is the signature of Ve appearance. 

Beam neutrino events in SK are selected based on 
timing information from the global positioning system 
(GPS) to reject cosmic-ray and atmospheric neutrino 
background. We start with the 112 fully-contained events 
in the 22.5 kiloton fiducial volume (FCFV) for K2K-I-HII. 
The definition of the electron signal candidates is given 
in Ref . |15| . The signal selection begins by identifying an 
electron ring candidate based on the shape of Cherenkov 
ring pattern and opening angle. Next, signal candidates 
are required to have electron-equivalent energy (visible 
energy) above 100 MeV to remove misidentified charged 
pions via v^-CC interactions and electrons from muon 
decays in which muons are below Cherenkov threshold. 
Finally, we require that no electron is followed by a muon 
decay within a 30 /isec time window. 

Near the limit of ^13 at the Aml^.^ region 

(sin2 26'^e ^ 0.05, Am^ ~ 2.8 x 10^^ eV^), we obtain 
an oscillation probability for u,, appearance in K2K of 
~ 10~^ , which gives an expectation of a few sig- 
nal events in our entire sample. Because this expected 
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appearance signal is so small, precise understanding of 
background and its reduction are crucial. One of the 
background sources is beam u^. This is responsible for 
13% of background events after the selection criteria de- 
scribed above. The remaining 87% originate from beam 
v^; this includes v^-CC and neutral current (NC) inter- 
actions of Vfj^ and i^^, where Vt appears as a result of 
oscillation. This f/^-originated background is dominated 
by tt" produced via NC interaction. A single 7r° decay- 
ing into two gamma rays can be classified by the stan- 
dard atmospheric neutrino reconstruction procedurei3| 
a single-ring showering event, when one gamma ray is 
not reconstructed, due to highly asymmetric energies or 
a small opening angle between the two gamma rays. 

For the current analysis, we enhance the tt'^ recon- 
struction capability with a newly introduced algorithm 
for electron candidate events. A second gamma-ray ring 
candidate is reconstructed by comparison of the observed 
charge and expected light patterns calculated under the 
assumption of two showering rings. The direction of the 
ring reconstructed by the standard procedure, visible en- 
ergy and vertex information are used as input. We use 
the invariant mass of tt*^ 77 {Minv) from the momenta 
of two potential gamma-ray rings. Among single elec- 
tron candidates, the 7r° background events have a value 
of Minv close to the 7r° mass. In contrast, for a Ve signal 
event, the reconstructed energy of the second fake ring 
tends to be small after sharing the visible energy with 
the two-ring configuration. Therefore Minv has a depen- 
dence on the energy of the electron candidate ring(i?e). 

For the final selection of candidate events, we exclude 
the region Minv > 100 MeV/c^. Then a tighter cut is 
applied to the events with Ee below 400 MeV, where the 
resolution of tt" mass is poor and the Ve signal shows 
small Minv The resultant Minv vs. Ee plot is shown 
in Fig. ^ The signal region is optimized to maxi- 
mize statistical significance assuming (sin^ 20^ej Am^) 
= (0.05,2.8 X lO-^eV^). With this requirement, 70% of 
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FIG. 1: Distribution of candidate events in the Mi„v-Ee 
space. The selection boundary is indicated by a line with 
a bend for (a) signal MC events and (b) z/^-originated back- 
ground MC events. Remaining events before the tt" rejec- 
tion in data (circles) are also shown for K2K-I and K2K-II 
combined. Signal MC distributions are obtained under the 
assumption of sin^ 26 = 0.05 and Am^ = 2.8 x 10"^ eV^. 



f^-originated background in electron candidate events is 
suppressed, and the overall efficiency for selection of Ve 
via CC interactions for K2K-I and K2K-II is 47% and 
51%, respectively. This selection efficiency of K2K-II is 
a little larger than K2K-I, however, the background 
expectation is also larger for K2K-II. The reduction of 
the data and background expectations according to the 
Monte Carlo (MC) simulation is summarized in Tabled 
After all selection cuts, we obtain one signal candidate in 
the data; four of five remaining before the final selection 
are rejected as 7r°-like. (The single signal candidate of 
K2K-I reported in Ref. [l5j is rejected by the new cuts.) 
This rejection capability is as expected, and a signal can- 
didate is consistent with background expectation. A vi- 
sual examination reveals the surviving candidate to be 
an event with more than two rings. It is consistent with 
a multihadron production event; these make up 14% of 
the background according to the MC simulation. 

Background expectation. — In our MC simulation, neu- 
trino interactions with oxygen nuclei are simulated as in 
Ref.p^. The interaction models used in this analysis 
are the same as the v^^ Vx oscillation analysis ex- 
cept that we set the cross section for CC coherent pion 
production to zero, based on Ref. [l7j . 

The energy spectrum and normalization are derived 
from measurements at the near detectors. The far/near 
V spectrum ratio used to extrapolate near detector mea- 
surements to SK is calculated using a beam MC simula- 
tion. In the simulation, we employ the tt"*" production 
cross-section measured by the HARP experiment pLSf. 
This spectrum ratio is validated by in situ measurements 
of pious from the aluminum target [lflj|. With the full 
data set the resulting flux is consistent with the previous 
result [^. The normalization of the events is deter- 
mined by the IKT data. The expected total number of 
FCFV events at SK without oscillation for K2K-I and 
K2K-II is 81.1 and 77.4, respectively. 

The number of background events is estimated to be 
0.8 (0.9) for K2K-I (K2K-II); 0.6 (0.7) event is originated 
from and the remaining 0.2 (0.2) is beam v^. The 



TABLE I: Reduction of events for search at SK. 

The first column lists each selection requirement. The others 
give, for each selection, the number of observed events, v^j, 
background in the case of no oscillation, and beam v^, for 
K2K-I and K2K-II, respectively. 



K2K-I K2K-II 



Data Vp, Beam-i^e Data Beam-i^e 



FCFV 


55 


81.1 


0.8 


57 


77.4 


0.9 


single ring 


33 


50.1 


0.5 


34 


49.4 


0.5 


shower like 


3 


2.7 


0.4 


5 


3.2 


0.4 


visible energy cut 


2 


2.5 


0.4 


5 


2.9 


0.4 


no decay electron 


1 


1.9 


0.3 


4 


2.2 


0.4 


non-7r° like 





0.6 


0.2 


1 


0.7 


0.2 



4 



background events are dominated by NC interactions and 
are consequently nearly insensitive to oscillated to Vr] 
Vt oscillation affects the background estimation by ^ 5% 
at the oscillation parameters of the K2K result |g| . 

Systematic uncertainties. — The various contributions 
to the systematic uncertainty in j/^-originated back- 
ground are summarized in Table |n| Most of the er- 
ror sources are related to tt" background in NC inter- 
actions. The fraction of NC single 7r° (NCItt") produc- 
tion in the total of background events in SK is about 
70%; the IKT detector data were used for understanding 
its cross section in water. The IKT measurement of the 
cross section ratio of NCItt" production to CC interaction 
(NCIttO/CC) gives 0.064±0.001(stat.)±0.007(sys.) [13, 
and our MC simulation predicts the ratio to be 0.065. 
Taking into account this measurement and the difference 
of the central values, we apply an uncertainty of 12% to 
the NCl7r°/CC ratio. For the composition of NC back- 
ground other than NCItt^, an uncertainty of 20% in the 
ratio of NC cross section to CC is quoted as in Ref . [2l| . 
For coherent NC 7r° production, we account for the dif- 
ference between the zero coherent pion production case 
and the model of Rein and Sehgal jl^]. Varying these 
cross sections by their errors, the contributions to the ex- 
pected number of events are estimated. In addition, final 
state interactions of nucleons and mesons inside of nu- 
clear matter could substantially alter 7r° momenta. This 
effect is estimated to be 8% on the i/^-originated back- 
ground, considering the differences in the tt'^ spectrum 
between IKT data and MC events. 

Systematic uncertainty related to the tt^ rejection has 
been evaluated based on the measurement of atmospheric 
neutrinos in SK. For this purpose we select two data 
samples which include an unreconstructed gamma ray 
from tt" decays in the final state. One of these is a sam- 
ple of single electron events; the other is a sample of 
events with one muon-like and one electron-like ring re- 
constructed. Using both samples, an uncertainty of 19% 
is estimated in the Minv cut in order to account for the 
effect of possible reconstruction biases which affect the 
Minv distribution. Water properties also affect the 7r° 
reconstruction and associated errors are estimated using 
cosmic ray muons. The light attenuation has been mea- 
sured during the entire period of operation of the SK 
detector. The absorption and scattering lengths were 
varied within their ranges accordingly, ±20%(±15%) for 
SK-I (SK-II), which results in systematic effects of ±11% 
(±6%) for K2K-I (K2K-II). 

For the heam-Ve contamination, the expectation is de- 
rived from the i^e/f/j flux ratio at SK with the beam 
MC simulation and flux extrapolation from the near 
detectors to SK. The beam MC expectation has been ver- 
ified by measurements of Vel^^i interaction ratio at the 
FGD complex. The measurement by SciBarJ3|j cover- 
ing the electron energy above 0.5 GeV, gives this ratio 
as 1.6 ± O.Z{stat.) ± Q.2{syst.)%. The ratio measured by 



TABLE II: Systematic uncertainties [%] in the expectation of 
!^^-originated background. When estimating the total uncer- 
tainty, the correlations between the neutrino fluxes and the 
cross sections arc taken into account. 







KOK TT 




fla") fib) 




NCl7r"/CC ratio 


±8 If? 


+ y 

— 7 


NC/CC ratio (non-NCl7r°) 


±4 ±3 


±3 


tt" energy spectrum 


±8 


±8 


coherent vr'^ model 


+3 
-10 


+3 
-10 


7r° mass cut 


+ 19 
-18 


+ 19 
-17 


water properties 


±11 


±6 


neutrino flux at SK 


1?? ±6 


±6 


non-QE/QE ratio 


+2 +1 
-2 -1 


+1 
-1 


detector efficiency 


±6 


±6 


single electron selection 


+5 
-7 


+7 
-8 


total 


+33 +37 
-32 -26 


+3y 

-24 



for the electron energy above 1 GeV is consis- 
tent with that of SciBar, and both show agreement with 
the beam MC prediction of 1.3%. The systematic un- 
certainty in the number of beam-t^e background is domi- 
nated by our understanding of pion and kaon production 
in proton collisions on the aluminum target; v,, contami- 
nation via muon and kaon decays is taken into account, 
where muons are emitted by positive pion decays. Uncer- 
tainties in the beam-z^e background of 14% and 16% are 
estimated for pion and kaon production, respectively. In- 
corporating other possible contributions such as the sin- 
gle electron selection, we quote an uncertainty of 
in total. 

For the number of appearance signal events, the to- 
tal estimated uncertainty is 15% with dominant sources 
coming from the energy spectrum and the cross section 
ratio of CC interactions other than CC-QE (non-QE) to 
CC-QE measured by the near detectors. 

Limits on mixing parameters. — To set an excluded pa- 
rameter region in a 2-flavor oscillation model, we adopt 
a confidence interval construction with the Poisson dis- 
tribution using the expected number of events and the 
observed number; these values are l.T^o ^ (in the case of 
no oscillation) and 1, respectively. The expected num- 
ber of events Afobs can be represented by the sum of two 
background components and an appearance signal as 

Nobs =S^^(Am^g, Sin^ 26'^e) ± Bbeam-u, 

±5^^(Am2^,sin2 26'^e), 

where B^^ is the number of electron candidate events in- 
duced by or oscillated to v^, and Bbeam-i^^ is that of 
beam-^'e. The number of — > u,, appearance signal, S^^ , 
and Bu^ depend on the probability of —f oscillation. 
For the contribution of oscillation to By^ , we 

assume the best-fit parameters of the K2K i/^ disappear- 
ance analysis 6] . In the calculation of the upper limit 
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FIG. 2: The upper bound on —> oscillation parameters 
at 90% and 99% C.L. The sensitivities of the K2K experiment 
for each C.L. are also indicated with dashed lines. 



on sin 20^e, the effects of systematic uncertainties are 
incorporated into the probability densities and the uni- 
fied ordering prescription of Feldman and Cousins (i^l is 
apphed. 

Figure 121 shows the upper bound on the oscillation pa- 
rameters for two flavor mixing, at the 90% and 99% con- 
fidence level (C.L.). Neutrino oscillation from i'^ to I'e 
is excluded at 90% C.L. in sin^ 26*^6 > 0.13 at Am^^ = 
2.8 X lO-^eV^ 

Reactor i7f, experiments provide complementary results 
to a search for non-zero ^13 with a z/^ beam; currently 
the only result for i^e appearance mode is provided by 
K2K. We note our resulting upper limit at 90% C.L. 
is sin^ 26*13 = 0.26 at Amj^ = 2.8 x 10"^eV^ assuming 
sin^ 26 — ^ sin^ 2^13 and Am^g ^ ^"^13- In the same 
Am? region, the most stringent limit from reactor exper- 
iments is sin^ 26*13 =0.1 by the CHOOZ experiment ^3 
and a weaker limit of 0.16 is reported by the Palo- Verde 
experiment [25l |. showing agreement with our result. 

In summary, the K2K experiment finished taking data 
in November 2004. Starting from June 1999, we accumu- 
lated data which corresponds to 9.2 x 10^^ POT. Com- 
pared to the previous search [l5j|. we improved both the 
statistics and the rejection of tt" backgrounds. As a re- 
sult, we find no evidence for neutrino oscillations in the 
Ve appearance mode. A single electron candidate is con- 
sistent with background expectation. We set an upper 



limit of sin^ 26'^e < 0.13 at 90% C.L. at the best-fit pa- 
rameters of the i^fi disappearance analysis. 
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